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We study the properties of two dimensional topological spin hall insulators which arise through 
spontaneous breakdown of spin symmetry in systems that are spin rotation invariant. Such a 
phase breaks spin rotation but not time reversal symmetry and has a vector order parameter. 
Skyrmion configurations in this vector order parameter are shown to have electric charge that is 
twice the electron charge. When the spin Hall order is destroyed by condensation of skyrmions 
superconductivity results. This may happen either through doping or at fixed filling by tuning 
interactions to close the skyrmion gap. In the latter case the superconductor- spin Hall insulator 
quantum phase transition can be second order even though the two phases break distinct symmetries. 
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Several recent theoreticalii^i^ and experimental^ pa- 
pers have discussed a phenomenon known as the quan- 
tum spin Hall effect in two dimensional insulators. In 
such an insulator an applied electric field leads to a quan- 
tized spin current in the transverse direction. Clearly 
the effect cannot occur if SU (2) spin rotation symmetry 
is preserved in the low energy theory of the insulator. 
Initial discussions of the spin Hall effect focused on situ- 
ations where the microscopic Hamiltonian has spin-orbit 
interactions which are such that only the component 
of the spin is conserved. Alternately one can contem- 
plate phases of matter where full SU{2) spin rotation 
symmetry is present at the microscopic level but is spon- 
taneously broken in the ground state in a manner that 
enables a spin Hall effect. Such phases were discussed 
recently in Ref. H. 

In this paper we study various aspects of the sec- 
ond kind of spin Hall insulator where the spin symme- 
try is spontaneously broken. The broken spin symmetry 
is characterized by a vector order parameter which we 
will denote N . Consequently in two space dimensions 
skyrmion topological defects are allowed. We show that 
the quantized spin Hall effect leads (with some restric- 
tions discussed below) to a quantized electric charge 2e 
on the skyrmion. Here e is the electron charge. This 
has some remarkable consequences. First so long as 
the skyrmion energy is much smaller than the energy to 
excite individual pairs of electrons, the skyrmion num- 
ber will be conserved due to charge conservation. Thus 
space-time hedgehog configurations of N are forbidden^. 
Destruction of the spin Hall order by condensation of 
skyrmions then very simply leads to a gapped s-wave su- 
perconductor. This may be done either by doping into 
the spin Hall insulator or by closing the skyrmion gap at 
fixed density by tuning interactions. In the latter case 
the resulting superconductor - spin Hall insulator phase 
transition can be second order despite their rather differ- 
ent broken symmetries. This provides a new example of 
a Landau-forbidden deconfined quantum critical poinl^. 

These results are best illustrated in the context of mod- 
els of of spinful electrons Cra hopping on a honeycomb 



lattice at half-filling: 
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where r, r' are nearest neighbor sites. Hint contains var- 
ious short-ranged interactions that preserve spin rota- 
tion and time reversal symmetries. Specific interactions 
that could stabilize the spin Hall insulating phase are 
discussed in Ref. |^. In the absence of Hint, the band- 
structure consists of two distinct Fermi points which can 
be chosen as K = (±47r/3>/3, 0) where the lattice con- 
stant is chosen as unity. Expanding the microscopic elec- 
tron annihilation operator Cr about these nodes in terms 
of continuum field operators Ciaa the low energy Dirac 
Hamiltonian is obtained as 
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Here i = L/R is a sublatticc index, a = +/— the node 
index, and a the spin. The corresponding Pauli matrices 
arc denoted x, fi and cr respectively. Further the trans- 
formation c+ = T^ip-f- enables writing the the Hamilto- 
nian in a manifestly symmetric form 
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The corresponding action (in real space) is 



5' ^ 



(3) 



(4) 



r being imaginary time, ^ 
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Here ^ = T,x,y with 

—iijj'^T^ and 7^,7^,7*' = jT'^ , —t^ . We have rescaled 
space-time to set the Dirac velocity to I. 

Consider a phase where the interactions Hint leads to 
a non-zero expectation value for the operator tpcrip, say 
along the z direction in the spin-space: 
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Within a mean field description this leads to a mass for 
the low energy Dirac fermions which is opposite for "f and 
i spin: 
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where the mass m = XN^ with A determined by the in- 
teractions that lead to this order. At the mean field level 
this is identical to the model discussed in Rcf. U for the 
quantum spin hall effect. It breaks spin rotation but not 
time reversal symmetry. In a sample with a boundary, an 
electric potential difference V applied between the right 
and left edge leads to a pair of propagating edge states 
with opposite spinii^ii^. Each spin species carries a non- 
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vanishing charge current of magnitude j^V transverse to 
the electric field in opposite directions. This corresponds 
to a spin-hall conductivity a^y = ^ (^7^) = e/27r. 

It is useful to characterize the spin hall effect in terms 
of the response to two different external gauge fields A° 
and A'^ which couple to spin and charge currents respec- 
tively. Consider the mean field action in the presence of 
these gauge fields 



(7) 

in the imits e = h = 1. When the fermions are integrated 
out, a non- vanishing transverse spin- hall conductivity of 
magnitude l/27r implies that the low energy effective ac- 
tion for the gauge fields is given by 
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Consider now the effect of fluctuations. These may be 
usefully discussed by considering the Dirac action in the 
presence of a fluctuating unit vector field N describing 
the orientation of the spin Hall order parameter: 
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We may now integrate out the fermions. In the limit of 
large mass m the result is a non-linear sigma model: 



5"= / d^x-^ [d^,N 
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with 'stiffness' - ~ |m|. The ellipses represent higher or- 
der terms in the 1/m expansion. Clearly in the ordered 
phase there will be two gapless linear dispersing Gold- 
stone modes associated with the broken spin symmetry. 

The other class of excitations associated with the vec- 
tor order parameter N are skyrmion configurations. The 
quantum numbers of a skyrmion whose size is much 
bigger than the length scale 1/m can be conveniently 
discussed within an adiabatic approach where such a 
skyrmion is slowly built up from the ground state. We 



now show that such a 'fat' skyrmion has electric charge 
2. 

To see this most simply, lets consider the Dirac action 
Eqn. [5] in the presence of a static background configura- 
tion of the N field that corresponds to a single skyrmion. 
An example is the configuration 

N{r,6) = {sin{a{r))cos{9), sin{a{r))sin{9)^cos{a{r))) 

(11) 

with the boundary conditions a{r = 0) = and a{r — s- 
00) = TT. Here (r, 9) are polar coordinates for two dimen- 
sional space. This field configuration corresponds to a 
skyrmion with Pontryagin index one. Rotate the N vec- 
tor field to the z direction by a unitary transformation 

U at all points in space: W {^a.N^ U ^ . Further, one 

may define tl.' = Uip' so that Eqn. [5] becomes 



d^x Vy (7^ (-ia^ + B,,) + imcr^) V' 
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where i?^ = —iU'^d^U . One readily finds that in the 
far field limit B{r ^ 00, 9) where eg is the unit 

vector along 9 direction. Therefore a skyrmion with unit 
Pontryagin index induces a spin gauge field A'^ — 2eo/r. 
The total flux of this gauge fleld is 47r which remains in- 
'Vvariant under smooth deformations of the N field (which 
does not change the skyrmion number). 

Let us gradually build in a skyrmion configuration by 
starting from the ground state. By the above argument, 
this threads in a spin gauge flux of 47r for A'' . Due to the 
spin-hall effect of the mean field state described above, 
this would result in a flow of electric current of magni- 
tude = in the radial direction. The total charge 
transferred from the center of the sample to the boundary 
during the process equals = rdO where the con- 
tour C is a circle located near infinity. This then leads to 
a charge = 2e associated with the skyrmion. Thus a 
skyrmion with Pontryagin index one carries an electrical 
charge 2e. 

This physical argument is confirmed by a more sophis- 
ticated analysis which carefully integrates out the Dirac 
fermion fields in the presence of an external charge gauge 
field. Specifically we consider 

S[A''] = J d^x^(7^ {-id^, + A^) ^/j + im'^a.N)ip (13) 

Integrating out the fermions and using the large mass 
expansion yields^i^ 



5" = 
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where Jj = ^e^^'^^N.d^N x dxN is the topological cur- 
rent whose time component equals the skyrmion density 
associated with the order parameter fieldiS. It follows 
that the skyrmion carries charge which is given by 



skyrmion 



= 2e / dxdy—e^^^N.d^N x dxN 
8n 
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Thus one again reaches the conclusion that a skyrmion 
with Pontryagin index 1 carries an electric charge two. 

It is important to emphasize that both the adiabatic 
argument or the field theoretic derivation above is really 
valid only for fat skyrmions with size much bigger than 
1/m. It is precisely such fat skyrmions that are impor- 
tant for our considerations below. 

The electric charge of the skyrmions has profound im- 
plications for quantum phase transitions out of the spin 
hall insulator. First if the energy of a fat skyrmion is 
much smaller than the energy of individual pairs of elec- 
trons, then the skyrmion number is conserved. Thus 
hedgehog configurations of the N vector in spacetime 
(which correspond to events that change the Pontrya- 
gin index) are prohibited. Further in the absence of any 
topological terms in the effective action for the N field, 
the skyrmions will be bosonsi^. Consider therefore the 
result of condensing the skyrmions by tuning the inter- 
actions in Hint while keeping the electron density fixed. 
Due to the charge on the skyrmions the result will be a 
superconductor! This superconductor will have a single 
electron gap, and fiux quantization in units of Indeed 
it is an ordinary s-wave superconductor. 

Remarkably the quantum phase transition between the 
spin hall insulator and the superconductor can be second 
order. This is despite the rather different broken sym- 
metries in the two phases. This is because the skyrmion 
condensation transition can simply be understood as the 
phase transition of an 0(3) vector model in the absence 
of hedgehogs. This is conveniently described in the CP^ 
formulation of the N field by writing N — z^az field with 
z a two component complex spinon field: 



(16) 



Here z — (zi,Z2) are coupled to an internal gauge field 
Ou = —iz'^dnZ. The skyrmion number Q is simply given 



by 
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Since the skyrmion number is proportional to the flux of 
the CP^ gauge field, of Eqn. [16] must be regarded 
as non-compact, and the phase transition is described by 
that in the non-compact CP^ modelii {NCCP^). 

Exactly the same field theory describes the deconfined 
quantum critical point separating Neel and Valence Bond 
Solid (VBS) phases of insulating square lattice quantum 
antiferromagnetsi. However there is an important differ- 
ence: for the Neel- VBS transition quadrupled hedgehogs 
are actually allowed but have been argued to be irrelevant 
at the critical fixed point of the NCCP^ model. They are 
however relevant at the paramagnetic 'free photon' fixed 
point. This leads to the appearance of two diverging 
length scales at the Neel- VBS transition. In contrast in 
the present problem the skyrmion number conservation 



is an exact symmetry at energy scales below the single 
electron gap which is large and finite in both phases and 
at the transition. Thus a suitable theory that captures 
the low energy physics of both phases and their transition 
should simply ignore all hedgehogs. Consequently there 
is only a single divergent length/energy scale at the de- 
confined critical point separating the spin hall insulator 
from the superconductor. 

The competition between the spin Hall order and 
the superconducting order can also be fruitfully de- 
scribed in a different manner which sheds further light 
on the role played by the quantum nature of the topo- 
logical defects near the phase transition. We com- 
bine the three component spin Hall order parameter 
N with the two component superconducting order pa- 
rameter Ipse into a five component 'superspin' vector 
(pa = {N'^ , , , Reipsc, IiTT-i'sc)- and impose a unit 
length constraint on (j)a, i.e X)a=i 'f'a^l'a = 1- What is the 
structure of a 'sigma model' that describes the dynamics 
of 0a? This question may again be answered by coupling 
the underlying Dirac theory to this superspin vector and 
integrating out the fermions. Specifically we consider the 
Dirac action 
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where the last term indicates the coupling to the super- 
conducting order parameter V'sc {ip'^ t"^ a"^ fi^ ip} ■ It is 
potent to make a unitary transformation on the fermions 
in the tp basis as V' = e^'^'^^^^ip and expand -0 in its Ma- 
jorana components ip = r]i + ir]2- In the rj basis the 
kinetic part of the Dirac action has manifest 5*0(8) sym- 
metry under rotations of all 8 components (2 spin, 2 node 
and 2 Majorona indices). The spin- hall order parame- 
ter is given as {rj'^T^a^r], rj'^r^a^ p^r], rj'^T^a^T]} while the 
s-wave SC has the form ^rj'^ t'^ fx"^ -q , tj'^ t'^ p,"^ r]} . 
Here p are Pauli matrices acting on the Majorana in- 
dex. These five components together transform as a vec- 
tor under an 5*0(5) subgroup of 50(8). This algebraic 
structure largely determines important features of the 
effective sigma model that describes the fluctuations of 
the (/)a when the fermions are integrated out. Within the 
large mass expansion the techniques of Ref. d lead to an 
50(5) action with a Wess-Zumino-Witten(WZW) term 
at level- 1. 
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The WZW term F is defined as follows. The field (pa 
defines a map from spacetime 5^ to 5^. F is the ratio of 
the volume in 5** traced out by pa to the total volume 
of 5^. Specifically let (p{x, u) be any smooth extension of 
p{x) such that 0(x, 0) = (1, 0, 0, 0, 0) and (p{x, 1) = (p{x). 
Then 



^abcde 

Area(5'4) ^ 



du / d^ xcpadxpbdypcdTpdduPe (20) 
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The presence of this 'topological term' crucially alters the 
physics from naive Landau-like descriptions of the com- 
petition between the two different orders. Parenthetically 
we note that the action above must be supplemented with 
some anisotropy between the N and ipse which breaks 
the 5*0(5) symmetry down to 50(3) x [/(I). The 50(3) 
corresponds to spin rotations and the U{1) to charge con- 
servation symmetry. Exactly the same sigma model with 
the WZW term also describes^^^ the Neel-VBS transi- 
tion of insulating quantum magnets in agreement with 
our previous identification. 

The physical implications of the WZW term are 
brought out by the following instructive calculation. 
Consider a hc/2e vortex defect in the superconducting 
state. In such a configuration the 5-component unit vec- 
tor will point along the 4, 5 directions far from the vortex 
core but will point along the 1, 2, 3 direction in the core. 
Specifically consider a configuration (we use polar coor- 
dinates (r, 6) for two dimensional space and r for time) 

4>{r,0,T,u) = {sin{a{r,u))N(T,u), 

cos{a(r, u))cos{9), cos{a{r, u))sin{dl^'^l) 

with N'^ = 1. Choosing a{r = 0,u) ~ 7r/2, a(r = 
oOjM) = Vu 7^ and a{r,0) ~ 7r/2, this configuration 
indeed describes a static superconducting vortex where (fi 
points along the N direction in the core. We have allowed 
the unit vector N in the core to have time dependence. 
The integrals defining T are readily evaluated for this 
configuration and lead to the result 

r = — / drduN.drN X duN (22) 
47r J 

This is precisely the quantum Berry phase for a spin-1/2 



object. Thus we see that close to the transition the super- 
conducting vortices behave as spinous. Condensing the 
vortices destroys the superconducting order but at the 
same time leads to condensation of the spin Hall order. 
This is the mechanism for the Landau-forbidden transi- 
tion. The arguments of Ref. [3 now establish the equiva- 
lence of the sigma model to the NCCP^ field theory and 
the superconducting vortices are directly identified with 
the CP^ spinous z. 

Finally we consider the result of doping the spin Hall 
insulator focusing again on the situation where the en- 
ergy of a fat skyrmion is much smaller than the en- 
ergy gap of individual electrons. Then the doping will 
be accommodated by the introduction of a finite den- 
sity of skyrmions into the system. At very low dop- 
ing a skyrmion crystal will presumably be stabilized but 
with increasing doping a translation invariant skyrmion 
condensate is expected. This will be a superconducting 
phase with no broken spin rotation symmetry. 

The route to superconductivity discussed in this paper 
is rather different from two common theoretical mech- 
anisms - phonon or other boson mediated pairing of 
electrons, or resonating valence bond pairing originat- 
ing from superexchange in a Mott insulator. Rather it 
involves condensation of charged solitons of a different 
broken symmetry, namely the spin Hall order. Similar 
phenomena have previously been considered in field the- 
oretic contexts^, and is known as topological supercon- 
ductivity. Superconductivity arising from the quantum 
spin Hall state thus provides a nice realization of a certain 
kind of topological mechanism for superconductivity. 
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